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ABSTRACT 

The U.S. Marine Corps deploys Marine Air-Ground Task 
Forces (MAGTFs) by airlift or sealift to participate in 
numerous short-term exercises. These exercises are of such 
duration that resupply of the MAGTF by strategic airlift or 
sealift is not practical. Thus, only stocked spare parts 
are available for repairs during the exercise. 

A model is developed which provides the operational 
commander with a stockage policy for spare secondary 
reparables (e.g., tank engines, amtrack transmissions, etc.) 
that optimizes the probability of successful mission 
completion subject to weight or volume constraints imposed 
by the MAGTF 's mode of deployment. Optimization of this 
stockage policy is stochastically modeled using data from 
the Marine Corps Integrated Maintenance Management System 
data base and then solved as an integer program. 

The integer program is coded using the Generalized 
Algebraic Modeling System language and solved using the 
Zero/One Optimization Methods mixed integer program solver. 
Operational data for a Marine Amphibious Unit yields an 
integer program with 190 binary variables and 26 
constraints. A solution within 0.07% of optimality is 
obtained on an IBM 3033AP computer in 3.9 seconds and on a 
Zenith Z-248 personal computer in 176 seconds. 
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I. 



INTRODUCTION 



The combat structure of the United States Marine Corps 
( USMC ) centers around the concept of the Marine Air-Ground 
Task Force (MAGTF) which is a combined arms force consisting 
of ground, air, and combat service support forces. The size 
of a MAGTF and its specific composition are dependent upon 
the particular mission assigned and the capability of the 
opposing forces. However, all MAGTFs are organized to take 
maximum advantage of the combat potential inherent in a 
rapidly deployable, closely integrated air-ground team under 
the control of a single commander. A MAGTF consists of four 
basis elements: command element, ground combat element, 
aviation combat element and combat service support element. 
It is the mission of the combat service support element to 
provide all logistical support for the ground and aviation 
combat elements. [Ref. l:pp. 1-4]. 

In this study, a model is developed which provides the 
operational commander with a stockage policy for spare 
secondary reparables (e.g., tank engines, amtrack 
transmissions, etc.) that optimizes the probability of 
successful mission completion within the constraints imposed 
by the MAGTF' s mode of deployment (i.e., airlift or 
sealift). The problem of determining this optimal stockage 
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policy is stochastically modeled using operational data 
resident in the Marine Corps Integrated Maintenance 
Management System (MIMMS) data base and then solved as a 
mathematical program. The remainder of this chapter will 
include a description of the scenarios which are of 
particular interest to the U.S. Marine Corps, an explanation 
of terms, and a definition of the problem. Additionally, 
justification for the approach taken here as compared to 
those in standard inventory models is provided. 

A. PROBLEM STATEMENT 

The U.S. Marine Corps currently deploys MAGTFs by 
airlift or sealift to participate in numerous short-term 
exercises such as BRIGHT STAR in Egypt, GALLANT EAGLE in the 
southwestern United States, and COLD WINTER in Norway. 
These MAGTFs range in size from 2,000 man Marine Amphibious 
Units (MAUs) to 16,000 man Marine Amphibious Brigades 
(MABs), depending upon the particular exercise. These 
exercises are of such duration that resupply of the MAGTF’s 
ground combat element by strategic airlift or sealift from 
logistics bases in the continental United States is not 
practical. Furthermore, by doctrine [Ref. l:pp. 1-4], a 
MAGTF is task organized to maximize combat power; therefore, 
the number of mechanics, technicians, and other logistics 
personnel together with their test and repair equipment are 
kept to a minimum. Thus, the quantities of spare repair 
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parts stocked by the combat service support elements are all 
that are available to support the requirements of the ground 
combat element for the duration of the exercise. 

In addition to these operational exercises, the Marine 
Corps deploys Marine Amphibious Units (MAUs) to the 
Mediterranean Sea as part of the U.S. Sixth Fleet on a 
continuous, six month rotating basis. These Mediterranean 
MAUs remain aboard ship except when conducting exercises. 
Therefore, the inventory of spares maintained by their 
combat service support element, called a MAU Service Support 
Group (MSSG), is limited by the cargo capacity of the 
assigned amphibious ships. Even though resupply from the 
continental United States by air is possible, it is desired 
to reduce dependency on this costly method [Ref. 2:p. 1] . 
Thus, similar to the above exercise scenarios, the MSSG must 
stock sufficient spares, within the capacity constraints of 
the assigned shipping, so as to satisfy the anticipated 
demands generated by the ground combat element's 
participation in exercises during the course of the six 
month deployment. 

Presently, a problem of even greater concern to Marine 
Corps combat service support organizations than either of 
those previously mentioned is the Marine Corps ' s role in the 
rapidly evolving maritime strategy of the U.S. Navy. As the 
amphibious power projection arm of the Navy, the Marine 
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Corps, in conjunction with the Navy, has developed and 
recently implemented the concept known as the Maritime 
Prepositioning Force (MPF). According to doctrine [Ref. 3], 
a MPF consists of a squadron of four to five ships loaded 
with combat equipment and supplies. One such MPF squadron 
is presently located in each of the following three ocean 
areas: Eastern Atlantic, Indian Ocean, and Western Pacific. 
When a MPF operation is ordered, a fly-in echelon consisting 
of 16,500 Marines and sailors, comprising a Marine 
Amphibious Brigade ( MAB ) and a Navy Support Element, is 
airlifted to a benign airfield in the vicinity of the 
objective area for link-up with equipment and supplies 
aboard the MPF squadron ships. The Brigade can be combat- 
capable and ready to move to an objective in five days or 
less and can operate for thirty days, independent of any 
strategic airlift or sealift from the United States, using 
the supplies and spare parts inventories on the ships [Ref. 
4:p. 11]. It is this last requirement, to operate 
independently for thirty days using organic supplies and 
inventories, that poses the greatest challenge to Marine 
Corps combat service support planners. However, as with 
supporting short-term exercises and MAU deployments, the 
problem of supporting the MPF MAB simplifies to one of 
stocking adequate spare parts inventories, subject to the 
capacity constraints of the MPF squadron ships, to meet the 
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demand during the MAB ' s projected missions. The fact that 
such demand may vary unpredictably must be confronted. 

At the present time, to solve spare stockage problems 
similar to those described above. Marine Corps combat 
service support planners determine the number of secondary 
items (referred to in Marine Corps' terminology as secondarii 
re^arabl_es ) to stock for each principal end item in the 
particular MAGTF's Table of Equipment (T/E) by using either 
last-period demand method or a simple moving average of past 
or historical demands generated during similar exercises or 
Mediterranean deployments. In the situation where no 
historical demands are available, such as with the new MPF 
MABs, stockage levels are determined using estimated 
replacement rates in combination with the experience of 
maintenance and supply personnel. These techniques have 
proved inadequate [Ref. 2:p. 1] because they fail to 
consider such factors as changes to a MAGTF's Table of 
Equipment and variations in the operational schedule, both 
of which may occur from one exercise or deployment to the 
next. More importantly, they neglect to consider that the 
demands placed on the supply system are the direct result of 
equipment failures. 

To overcome these inadequacies, it is the objective of 
this study to develop a model which provides the operational 
commander at the MAU or MAB level with a stockage policy for 
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spare secondary reparables that optimizes the probability of 
successful mission completion within the constraints imposed 
by the mode of deployment. The proposed model is 
sufficiently general so as to be adaptable to any of the 
above scenarios and it can accommodate any additions or 
deletions to the MAGTF's Table of Equipment or variations 
in mission duration among any combination of principal end 
items. Additionally, the model incorporates a measure of 
effectiveness which provides the operational commander with 
an indicator of the overall impact of the stockage policy on 
successful mission completion. Furthermore, the model 
utilizes the operational maintenance data resident in the 
Marine Corps Integrated Maintenance Management System 
(MIMMS) to determine operationally based failure rates. 
Thus, by use of this model a nearly optimal stockage policy 
for secondary reparables at the retail level can be 
determined for a wide variety of inventory problems 
currently facing Marine Corps combat service support 
planners. This will result in a more effective utilization 
of the limited space and weight resources available for 
spares' inventories on shipping or aircraft assigned to 
deploying Marine Corps units. 
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B. INVENTORY MODELS AND SOLUTION APPROACHES 



The stockage problems confronting Marine Corps planners 
can all be viewed as variations of the single period, 
constrained, multi-item inventory problem. The simplest 
such single-period model is the classic newsboy problem 
[Ref. 5], alternatively referred to as the Christmas tree 
problem. The characteristic distinguishing the models of 
this thesis from other standard inventory models, such as 
the classical economic-order-quantity or economic-order- 
interval models [Ref. 6:pp. 159-289] is that only a single 
time period, usually of finite length, is relevant; 
therefore, only a single initial procurement or order is 
made [Ref. 6:p. 297]. The previously described scenarios 
all fall within the context of single-period models because 
the demand for secondary reparables is the direct result of 
equipment failures which occur at infrequent intervals 
during the course of the exercise. Resupply either does not 
occur or is kept to a minimum, and thus it can be ignored. 
Furthermore, it is Marine Corps policy to permit only 
functional unit replacement at the field or organizational 
levels of maintenance. Actual repair of secondary reparables 
at the intermediate maintenance level is very limited and 
for all practical purposes performed only at the depot 
level. Therefore, since the MAGTFs of concern in these 
scenarios are usually only authorized up to the intermediate 
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level of maintenance, the actual repair and return of 
secondary reparables to these retail levels of inventory is 
insignificant and can be ignored. This is particularly true 
when one considers that some of the exercise durations are 
of the same order of magnitude as the average repair time 
for some of the secondary reparables. Thus, each scenario 
consists of an exercise (or exercises) which is (are) of 
finite duration and during which repair and resupply are 
essentially nonexistent, so the combat service support 
planner must make a one-time decision of what to stock, not 
how much of each item to procure. He must do so within the 
constraints of the mode of deployment. This is analogous to 
the single period, constrained, multi-item inventory 
problem. 

There are several examples in the literature of single 
period inventory models applied to military logistics 
situations. They include the spare parts kit problem [Ref. 
7:pp. 281-295], the fly-away kit problem [Ref. 6:p. 324], 
and the submarine provisioning problem [Refs. 8:pp. 236-243, 
9:p. 407]. The problem with each of these military 
applications is that they assume that there is some monetary 
stockout cost incurred when a demand exists and no spares 
are available to meet the demand. This is realistic in the 
civilian sector where real monetary values can often be 
associated with each stockout. But, in most practical 
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military applications it is difficult, if not impossible, to 
determine a tangible, realistic, and objective stockout 
cost. This is particularly true for each of the previously 
described scenarios because there is no way to place a cost 
on the inability to complete a mission because of a 
deficiency of spares. Thus, it is difficult to apply the 
models existing in the literature to these real-world Marine 
Corps inventory problems. 

This problem is not insurmountable. In the civilian 
sector when an organization does not know its stockout costs 
or is not confident in its estimates of them, management 
will set service levels which indicate the ability to meet 
customer demands from stock [Ref. 10 :p. 149]. A service 
level is basically the complement of the probability of 
stockout and really represents a subjective decision by 
management to accept some degree of stockout or customer 
disservice [Ref. 10:p. 150]. The U.S. Navy actually uses a 
similar concept when establishing its Aviation Consolidated 
Allowance Lists (AVCAL) and Coordinated Shipboard Allowance 
Lists ( COSAL ) [Ref. ll:pp. 4-30-4-50]. The service level 
measure is less than satisfactory for military applications 
since it is subjectively established by combat service 
support planners. Service levels do not provide the 
operational commander in the field with a meaningful 
indicator of the impact of the recommended stockage policy 

9 
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on successful mission completion. This is because there is 
no direct relationship between the probability of stockout 
and the probability of successfully completing a given 
mission. A high service level does not equate to an optimal 
probability of mission completion. This is particularly 
true in the constrained, multi-item inventory situation 
where it might actually be better to accept a greater 
probability of stockout (lower service level) for a 
secondary reparable with a low failure rate in order to be 
able to stock more of another secondary reparable with a 
high failure rate. 

To overcome the disadvantages associated with the 
minimization of stockout costs (actually variable costs) and 
service level approaches, the proposed model introduces a 
new measure of effectiveness, the grobgbllgtg o£_ mls_slon 
completion given the specified stockage goll£y_, which will 
be referred to as the £robabi_lltg of mission cgmgle_tlon 
throughout the remainder of this paper. This probability 
extends the familiar concept of reliability to include a 
system and its inventory of spare parts [Ref. 8:p. 237], A 
system is defined to be the quantity of a principal end item 
required to be operational throughout the duration of a 
given mission. This end item must be a "readiness 
reportable pacing item" as delineated in Marine Corps 
Bulletin 3000 [Ref. 12:pp. 1-12 Enel (2)]. A readiness 
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reportable pacing item is a principal end item whose 

operation is essential to the ability of a MAGTF to complete 

its assigned mission. Simply stated, the probability of 

mission completion is the probability that k-out-of-n (k£n) 

readiness reportable pacing items of the same type remain 

operational for a specified mission with a given stockage 

policy. Black and Proschan [Ref. 7:p. 283] proposed a 

similar measure of effectiveness in their early work using 

marginal analysis to determine optimal spare parts kits at 

minimum cost. They called it pr obablllty qf_ adequacy or 

assurance of adequacy, which is defined as 

The probability that either the equipment does not fail 
due to a. . .malfunction or, if it does fail for this 
reason, a replacement is available to provide continued 
successful equipment operation. [Ref. 8:p. 237] 

However, this definition applies only to series systems in 

which all equipment must remain operational for a successful 

mission. The proposed probability of mission completion is 

more general since it can be extended to "k-out-of-n 

systems". This is required in the present study because 

the operational commander, not the military logistician, 

determines the quantities of a principal end item which he 

feels are necessary to complete the assigned mission. For 

example, a MAU commander may feel that he can accomplish a 

particular mission with just three out of five tanks which 

permits fewer tank secondary reparables to be stocked than 

if he required all five tanks to complete the same mission. 
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In this study the inventory problems in each of the 
previously described scenarios are modeled as single period, 
constrained, multi-item inventory problems in which the 
demands are the result of independent failures of secondary 
reparables. In the mathematical formulation of these 
models, the probability of mission completion is maximized, 
c as opposed to the standard approach of minimizing expected 

(variable) costs. This provides models which reflect the 
fact that in most military situations, and specifically in 
the three Marine Corps scenarios studied here, stockout or 
shortage costs cannot be determined. The inclusion and 
numbers of constraints are dependent on the mode of 
deployment. However, in all cases studied, they are linear 
and will be referred to as capacity constraints. Thus, each 
model will be viewed as a mathematical program in which the 
objective function to be maximized is a nonlinear, 
nondecreasing function of the decision variable representing 
the quantity of secondary reparables to stock. Inclusion of 
the linear capacity constraint or constraints result in an 
overall integer nonlinear programming problem. Since the 
objective function is itself a complicated function of the 
decision variable, the solution approach is to linearize the 
objective function and then reformulate the entire problem 
using binary variables. This reformulation converts the 
original integer nonlinear programming problem into a binary 
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linear programming problem which is then coded for computer 
processing using the Generalized Algebraic Modeling System 
(GAMS) [Refs. 13, 14, 15] and solved using the Zero/One 
Optimization Methods (ZOOM) [Ref. 16] mixed integer program 
solver . 

Since one of the primary purposes of this study is to 
provide the combat service support planners of the Second 
Force Service Support Group at Camp Lejeune, North Carolina 
with models solvable on a personal computer, it is desired 
to avoid dependency on commercial optimization packages not 
locally available. Therefore, the initial solution 
methodology utilized Langrangian relaxation to establish an 
upper bound on the optimal solution of the integer linear 
problem. Then a heuristic based on marginal analysis is 
used to improve upon the best feasible solution obtained 
during the establishment of this upper bound, and this 
heuristical ly improved solution is accepted as the final 
solution provided that its associated objective value was 
reasonably close to the upper bound. This methodology is 
similar to that implemented by DeWolfe [Ref. 17:pp. 28-32] 
to solve a Selective Reenlistment Bonus problem for the U.S. 
Marine Corps. The advantage of this approach is that it can 
be programmed in a language such as FORTRAN 77 thus avoiding 
dependency on a commercial optimizer. However, as the number 
of constraints increase the programming difficulty 
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associated with the search procedure used in establishing 
the best upper bound increases, and for any more than two 
constraints actually becomes prohibitive. So, to keep the 
solution approach as flexible and general as possible, it 
was decided to solve all integer linear problems with the 
GAMS/ZOOM optimization package ( commerically available in 
March 1988), thus allowing the solution of models with more 
than one constraint and avoiding restrictive programming 
considerations . 

The above solution approach overcomes the disadvantages 
associated with more standard solutions of the single period 
inventory model. For instance, it could be argued that 
since these inventory problems are all examples of decision 
making under risk, they could be solved using decision 
theory and payoff matrices to minimize expected costs. 
However, as was pointed out earlier, realistic monetary 
stockout costs are extremely difficult to determine in most 
military applications. Additionally, in the multiple item 
inventory problem, the payoff matrices become very large and 
the interrelationship among the payoffs associated with 
different combinations of items becomes difficult to 
evaluate and enumerate, so it is computationally more 
efficient to use a marginal analysis approach [Ref. 10:p. 
269]. Marginal analysis provides the basis for the 
techniques used in the spare parts kit problem [Ref. 
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7 : pp . 283-284] and in the submarine provisioning problem 
[Ref. 8:p. 237], but it has practical limitations when more 
than one constraint must be considered. Hadley and Whitin 
[Ref. 6 : pp 328-330] obtained exact integer solutions to 
their single constraint, fly-away kit problem using dynamic 
programming and Dreyfus and Law [Ref. 18:pp. 107-118] offer 
dynamic programming algorithms to solve single constraint, 
cargo-loading problems. However, they all point out that 
for problems involving more than two constraints even 
dynamic programming has computational limitations, 
particularly when both the number of constraints and 
decision variables are large [Refs. 6:p. 331, 18:pp. 107- 
118] . The solution approach utilized in this study avoids 
these shortcomings. 

C. THESIS OUTLINE 

This thesis develops and presents a method for 
determining spare stockage levels of secondary reparables 
associated with the readiness reportable principal end items 
of a given Marine Air-Ground Task Force. In Chapter II, the 
maintenance data base, data analysis, and problems 
encountered during this analysis are discussed. The 
objective functions are defined, developed, and justified in 
Chapter III. Additionally, each scenario is modeled and 
formulated as an integer nonlinear programming problem for 
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which the objective function is linearized so that the 
problem can be reformulated as a integer linear problem to 
facilitate its solution. Specific solution methodology is 
presented in Chapter IV. The initial approach using 
Lagrangian relaxation with a heuristic is briefly described 
before the more general procedure using GAMS/ZOOM is 
developed and implemented. In Chapter V results using 
actual data for a typical Marine Amphibious Unit are 
presented and verification of the model is conducted using 
different sets of simulated failure data. The chapter 
concludes with a discussion of computational experience for 
both the mainframe and personal computer. Chapter VI 
contains conclusions and recommendations. Finally, listings 
of both FORTRAN 77 and GAMS source codes are provided in the 
appendices . 
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II. DATA ISSUES AND ANALYSIS 



A primary impetus for this study was the desire to 
develop a reliability-based, spares stockage model for 
secondary reparables utilizing the maintenance data resident 
in the Marine Corps Integrated Maintenance Management System 
(MIMMS). Currently, spare stockage models within the Marine 
Corps are limited to those based on forecasts from past 
demands generated against the supply system. This 
forecasting approach neglects to consider that the source of 
the demands is the failure of equipment such as secondary 
reparables, and it is precisely this type of failure 
information that can be extracted from relevant MIMMS sub- 
files which permits the development of a reliability-based 
model . 

A. DESCRIPTION OF DATA 

MIMMS is an automated system which accumulates 
maintenance data on all serialized, principal end items in 
the Marine Corps inventory. It records and accumulates, in 
various sub-files, maintenance actions performed on a 
serialized end item during its service life at each of the 
three levels of maintenance (i.e., organizational, 
intermediate, and depot) [Ref. 19]. In particular, the 
MIMMS sub-file of relevance to this study is the Equipment 

17 



/ 



REPRODUCED AT GOVERN MTN I F CTFI'^I 



Repair Order History File. This sub-file maintains all 
equipment repair orders ( EROs ) that were completed on every 
serialized end item during the past eighteen months, and it 
is updated on a quarterly basis. An ERO is opened each time 
a serialized end item is inducted into the maintenance cycle 
and contains information such as status of the end item 
(i.e., combat deadlined, safety deadlined, mission degraded, 
etc.), the meter reading at the time of the current 
maintenance action (expressed in either hours, miles, rounds 
or days depending on the end item), and the defect code 
which describes the nature of the failure. By looking at 
only those EROs in the MIMMS ERO History File which were 
opened for critical, combat deadlining repairs (i.e., 
category code M), a new, local user file can be created. 
This user file is organized by principal end item and 
includes only those end items which are considered 
"readiness reportable" [Ref. 12:pp. Enel (1) 1-23]. 
Associated with each end item are all serial numbers for 
which a combat deadlining ERO exists; each serial number has 
a meter reading and defect code to reflect each ERO 
completed on that particular serial number during the 
eighteen month period covered by the ERO History File. 
Obviously, there may be multiple entries under the same 
serial number indicating that this end item suffered more 
than one combat deadlining failure during the reporting 
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period. There 


may 


be no entries 


if no failures were 


experienced 


by 


the 


particular 


serial number. Table 1 


provides a sample extract from 


this 


user-created data file 


and displays 


its 


contents for several 


tank serial numbers. 








TABLE 


1 






EXTRACT 


FROM LOCAL 


USER : 


DATA FILE 








METER 




DEFECT 


SERIAL # 


EOTC* 


READING 




CODE** 


502814 


H 




74 


D05 


(POWER TRAIN, CLUTCH) 


502828 


H 




51 


A34 


( ENGINE REPLACE ) 




H 




75 


J34 


(COOLING SYS REPLACE) 




H 




107 


C34 


(POWER RACK REPLACE) 




H 




155 


A3 4 


(ENGINE REPLACE) 




H 




163 


B34 


(TRANSMISSION REPLACE) 


502831 


H 




81 


A34 


(ENGINE REPLACE) 




H 




149 


B34 


(TRANSMISSION REPLACE) 


502835 


H 




62 


KOI 


(ELEC SYS, GENERATOR) 



* EOTC= Equipment Operating Time Code 

(H=Hours, M=Miles, R=Rounds, D=Days) 

** Defect Codes are alphanumeric codes and have been 
annotated here for display purposes only. 



The data in this consolidated file can be used to estimate 
failure rates for each secondary reparable in either of two 
ways and both methods will be described in the following 
section along with their advantages and disadvantages. 
It must be stated that both the format and method of update 
for many of the MIMMS files require improvement before a 
truly comprehensive analysis of the data is possible. A 
case in point is the ERO History File used in this study, 
and recommendations for its improvement are provided in the 
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following paragraph and in Chapter VI which minimize the 
impact on current MIMMS file management and the amount of 
effort required at the using unit input level. 

One of the basic problems with the ERO History File is 
that potentially useful data is being discarded during each 
quarterly update in order to minimize file size. What is 
proposed here is to maintain a local (i.e., at the Force 
“ Service Support Group level ) consolidated version of the ERO 

N' 

u. History File with a format similar to that described in the 

preceding paragraph. This file should be updated when its 
parent ERO History File is updated at the end of each 
quarter. Since this consolidated file contains only 
information extracted from combat deadlining EROs, which 
comprise a small fraction of the total EROs in the parent 
file, there is a greater capability to accumulate relevant 
data without the need to automatically discard the oldest 
three months of data. The basic algorithm for each 
quarterly update is: 

Step 1. Select a combat deadlining ERO (Category Code M) 
from the ERO History File if one exists. If none 
exist, STOP. 

Step 2. Determine if the serial number associated with 
this ERO is already resident on the consolidated 
file. If it is, continue with Step 3. If it is 
not, add the serial number and its associated 
meter reading and defect code to the consolidated 
file under the appropriate principal end item. 
Return to Step 1. 
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Step 3. Determine if the same defect code exists for this 
serial number. If it does, continue with Step 4. 
If it does not, add the meter reading and defect 
code to the list associated with this serial 
number. Return to Step 1. 

Step 4. Subtract the two meter readings to obtain a time 
between failure for this type of defect. Add 
this new data point to the list associated with 
this type of defect. Delete the record 

associated with the older meter reading (i.e., 
the smallest meter reading) and add the new meter 
reading and defect code to the end of the list 
associated with this serial number. Return to 
Step 1. 

The result is an updated version of the consolidated file 
which is retained at the local level for update during the 
next quarterly cycle. Each readiness reportable principal 
end item in the file is broken down by serial number which 
in turn contain records of the most recent meter readings 
and defect codes. Each end item also consists of lists 
containing times between failures for each type of defect. 
Thus this local consolidated file does not alter any 
existing MIMMS file definitions, nor does it necessitate the 
creation of any new system file definitions. More 
importantly, it provides the combat service support planner 
with a consolidation of all the pertinent information 
required to compute failure rates and track the life history 
of any secondary reparable associated with a readiness 
reportable end item. This is a significant improvement over 
the present method which is restricted to using the most 
recent eighteen months of failure data resident in the ERO 
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History File, and it is this restriction that limits the 
accumulation of failure times for certain types of secondary 
reparables (e.g., certain radio transmitters). 

B. DATA ANALYSIS 

Using the data resident in the consolidated file, the 
operating times between failures for each type of secondary 
if reparable were computed. Failure rates were estimated using 

two different approaches. In the first approach, an 
exponential lifetime distribution was assumed for each 
secondary reparable and the maximum likelihood estimator 
(M.L.E.) for the parameter / t - of an exponential distribution 
was used to provide an estimate of the failure rate. The 
second approach involves a Bayesian analysis of the times 
between failure to determine a posterior distribution for 
the failure rate of a secondary reparable. 

In deriving the M.L.E. for , it is assumed that the 
observed lifetimes of each secondary reparable i form a 
relevant random sample of size n^ from an exponential 
distribution with unknown parameter The likelihood 

function must also account for the exposure times during 
which no failures of type i occurred. If these terms are 
not considered, the computed failure rate will tend to be 
too large and excessive stockage levels would result. For 
the n^ observed failure times and the observed exposure time 



22 



['i during which no failures occurred, the likelihood 
function is: 




/.[ exp • ... • a i exp /mitin i • exp ' 



) n i 



exp 



f 

+ c 'i) 



>= i 



Maximizing the likelihood function with respect to 



yields the following estimator of the failure rate: 



rt; 



n ; 



'■i = 



V i- j. f ■ 

l_ji] • i 

y=i 



(2.1) 



where t,- represents the total exposure time for secondary 
reparable i. [Ref. 20:pp. 282-296] 

The Bayesian approach also assumes that the distribution 
of the lifetimes for each secondary reparable i is 
exponential with parameter / however, it further assumes 
that the exact value of this parameter is a realization of a 
random variable. Even though ^ is not precisely known, it 
has a prior distribution which is taken here to be a 
conjugate Gamma distribution with shape parameter 
a i(2i>0) and scale parameter // i > 0): 






A>-i ; «/ . Pi) = < 



PV 

r(« f ) 



A: 



1 1 exp 



A > 0 



(2.2) 



o o . 

The prior distribution parameters are estimated using a 
parametric empirical Bayes ( PEB ) approach for failure rates 
as outlined by Gaver and Lehoczky [Ref.21:pp. 220-224]. The 
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PEB approach uses the entire data set to compute estimates 

A A 

a j_ and which are substituted into the applicable 

formulas. To compute the prior parameters in PEB requires 
the maximization of the "marginal likelihood function" : 






> Pi > Hi ' Lj) — | | f 
m=l ( 



■ oo (/ . (. ) n im 



mrnn * 



n im 



exp /,m ^ m 






r(*i) 






S-'exp-^ 



im 



im d) • 



(2.3) 



where, 

m individual source of data for component i (e.g., 
tank serial #502826) 

M t total number of data sources for component i 

). im failure rate for component i of data source m 

hm observed exposure time of component i on data source 

m 

n lm number of failures of component i on data source m 

a / shape parameter of Gamma prior distribution of 

failure rate for component i 

Pi scale parameter of Gamma prior distribution of 
failure rate for component i. 

Integration results in the following marginal likelihood 
function for secondary reparable i: 



M; 



> Pj > Q-i > Li) — j~J 



m— 1 



r ( n im + «,) 
m ■" (“/) 



Pi 



im 



+Pi 



hm 



hm 



+ P 



\ n im 
i ' 



(2.4) 



which is the product of individual negative binomial 



distributions with parameters p= 



Pi 



l im + Pi 



and 



r = «/ 
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The numerical maximization of this function subject to known 
upper and lower bounds on each a and is discussed in 
Chapter IV. The important result of this maximization is an 

A A 

estimate for the shape ( aj_)and scale (fi±) parameters 
of the Gamma prior distribution that incorporates the actual 
observed operational data for each secondary reparable. 

Returning to the Bayesian analysis, the n^ observed 
times between failure and the observed exposure times i\ 
during which no failures occurred are all a function of 
A, and give rise to the likelihood function. 






.Tit 



= /.-exp 



- / r T / 



(2.5) 



where r i once again represents the total exposure time for 
secondary reparable i. We are now interested in the 
updated distribution of / ^ after the times between failure 
have been observed so as to provide a current estimate of 
future failure events. This is the posterior distribution 
and it is proportional to the product of the likelihood 
function (2.5) and the prior distribution (2.2): 



/(/•; i Lj) oc /(/,- I . 

Thus, 



/('’•/ U) = c^'exp-^0(A“-‘e.xp-^) 
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where is a constant factor which depends on the observed 
times between failure, but does not depend on X ^ . The value 
of this constant can be determined using normalization. The 
resulting posterior distribution is: 



A , A 

( T . 4 - P)) n i + a i A a 

A'-i I Li) = A'-i I data) = — — — + “* ” 1 exp” (T; ' + 



(2.6) 



rK- + «/) 



a Gamma distribution with shape parameter and scale 

A 

parameter t z - -f /? ; - . [Ref. 20 :pp. 257-280] 

The posterior distribution could also have been derived 
directly using the fact that the family of Gamma 
distributions serves as a conjugate family of prior 
distributions for samples from an exponential distribution. 
Thus, using the theorem for Gamma conjugate priors [Ref. 
20:pp. 271-272], the posterior distribution for ^ given 

the observed times between failure, is a Gamma distribution 
with shape parameter n^ + and scale parameter . 

This agrees with the result derived in the previous 
paragraph and indicates that the posterior distribution can 
be defined knowing only the number of observations and the 
sum of all the exposure times. This is an important result 
because it enables the posterior distribution to be updated 
very easily as new observations are obtained. 

Since the posterior distribution of A z - has been shown 
to be a Gamma distribution, its expected value is: 
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E[/., | data] = 



T i + Pi 

This expected value also happens to be a "Bayes estimator" 
of A when a "squared error loss function" is used [Ref. 
20:p. 227]. This estimator could also have been used to 
provide an estimate for each A j_. However, the primary 
purpose of the Bayesian analysis was not to provide a point 
estimate of A a, but rather to provide a distribution (i.e., 
posterior distribution) for each A • which could be used to 
make probabilistic predictions about demand functions 
involving A/. This step is detailed in the next chapter. 
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III. MODEL FORMULATION AND DESCRIPTION 



In this chapter, the two different objective functions 
considered in this study are defined, derived, and 
justified. Additionally, each of the three scenarios 
described in the introduction is formulated as an integer 
nonlinear program whose objective function is subsequently 
linearized to create an integer linear problem. 

A. MODEL DEVELOPMENT 

The following formulation is developed to model the 
determination of secondary reparable inventory levels for 
each of the three scenarios described in the introduction. 
With the exception of the number of constraints, the 
formulation is the same for each scenario. 



INDICES: 

/= 1,- ,iV 

DATA: 

w i 

W 

v i 

V 

“i 



secondary reparable or component 

unit shipping weight of secondary reparable i 
[lbs] 

weight capacity of aircraft or ship [lbs] 

unit volume of secondary reparable [cuft] 

volume capacity of aircraft or ship [cuft] 

maximum allowable number of spares for 
component i 
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FUNCTIONS: 



Pi(xi ) probability that the number of readiness 

reportable pacing items (specified by the 
operational commander) containing secondary 
reparable i will survive a mission given 
spares are stocked 

DECISION VARIABLE: 

x i Number of spares to stock for secondary 

reparable i 



FORMULATION: 

4 V 

max 

j=1 

N 

s.t. /WjXj < IV 

i = 1 

.V 

y < v 

i= 1 

o < X; < 10 ,.v 

X; integer 



(PI) 



The objective function is a mathematical expression of 
the probability that a quantity (specified by the 
operational commander) of each readiness reportable pacing 
item in a MAGTF's Table of Equipment will survive a mission 
of specified duration given a particular stockage policy for 
secondary reparables. This probability was referred to in 
the introduction as the probability of mission completion. 

29 
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In terms of this definition, successful mission completion 
hinges on the survival of the MAGTF's pacing items. Each 
pacing item is composed of secondary reparables. The 
failure of any one secondary reparable will cause the pacing 
item to be combat deadlined and unable to complete its 
mission. Thus, each pacing item is viewed as a system of 
secondary reparables in series. However, in a MAGTF's Table 
of Equipment, there is not just one of each type of pacing 
item, but rather a specific quantity of each type necessary 
to accomplish the assigned mission. This minimum quantity 
required will not only vary from one mission to another, but 
also among different MAGTF commanders when confronted with 
identical missions. Therefore, the different pacing items 
taken as a whole do not constitute a simple series system in 
which all components must survive the given mission. 

In defining the objective function, two different 
functional forms were derived and studied. Both rely on the 
assumptions that secondary reparables fail at constant rates 
and that the failures of individual secondary reparables 
are independent of one another. The independence assumption 
is certainly true for the secondary reparables associated 
with different types of pacing items; however, it may not 
necessarily be true for the secondary reparables of a 
specific pacing item. For example, the failure of M60A1 
tank engines is independent of the failure of LVTP7A1 
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amtrack engines, but the failure of amtrack final drives may 
not be independent of amtrack transmission failures. These 
assumptions together with the fact that secondary reparables 
from one type of pacing item are not interchangeable with 
those of another (e.g., a tank engine cannot be substituted 
for an amtrack engine) make it possible to consider only 
secondary reparables without having to specify the pacing 
items of which they are components. 

The functional form of the first objective function, in 
which the failure rate of each secondary reparable is 
assumed to be known, is as follows: 



pm = 



a: + X: — k: /s 

v h£i!A 



m 



/_ 

m—0 



ml 



exp 



A 

-Wi'i 



( 3 . 1 ) 



where , 

a i number of pacing items containing secondary reparable 
i 

ki minimum number of pacing items containing secondary 
reparable i required to be operational throughout 
given mission (/:, <«,-) 

A 

A/ M.L.E. of failure rate for secondary reparable i 

r i mission duration during which secondary reparable i 
must operate. 

This function is a direct consequence of the constant 
failure rate assumption for each secondary reparable. 
Underlying this constant failure rate is an inherent 
exponential lifetime distribution which implies that the 



31 



REPRODUCED AT GOVERNMF N I 



number of failures defines a Poisson (arrival) process with 
exponential interarrival times. The number of failures of 
secondary reparable i to occur in any fixed interval of 
time, such as a mission of duration tj_, will have a Poisson 
distribution with mean However, since there are a 
total of aj_ pacing items containing secondary reparable i 
and operating during this mission, the total number of 
failures will be the sum of independent, Poisson random 
variables (each with a mean jt ^ ) which is itself a Poisson 
random variable with a mean a,/.^ . This is really only 
an approximation since not all a_£ pacing items may be 
required during the mission. In which case the overall 
failure rate would be smaller implying that we are probably 
overestimating the number of spares required. On the other 
hand, it can be argued that with fewer than a^ items the 
overall failure rate of the operational items may increase 
since the same mission must be accomplished with fewer 
assets. In any case, this can be viewed as a mild, worst 
case approximation. 

Therefore, equation (3.1) represents the probability 
that there are less than or equal to some specified number 
of failures of secondary reparable i during a mission of 
duration tj_. This specified number is the sum of the number 
of spares, Xj_, and the difference between the number of 
pacing items containing secondary reparable i and the 
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required minimum number of pacing items containing secondary 
reparable i (k^) . As an example, consider the situation 
where the Table of Equipment of some MAGTF consists of five 
tanks each containing one engine (a^=5) and assume that the 
MAGTF commander decides that to accomplish the assigned 
mission he must have a minimum of four tanks operational at 
all times (k^ = 4). Then, assuming one spare tank engine 
(x^=l) is stocked, there can be no more than two engine 
failures (a^+x^-k^= 5+l-4=2) to adhere to the commander's 
guidance. It can be shown that this probability is equal to 
the probability that the quantity of each type of pacing 
item (or equivalently all its component secondary 
reparables) remains above a specified level (k^), because 
the following relation must always hold true: 

{number of failures of secondary reparable i} + 

{number of secondary reparables i operational and in 
inventory} = 

{total number of secondary reparables i, ( a ± + x ± ) } 

Thus, 

P( {number of secondary reparables i operational and 
in inventory at t}>k^) = 

P( (aj+Xj } - {number of failures of secondary reparable 
i in t } >_k^ ) = 

P( {number of failures of secondary reparable i in t} 
: <a i +x i -k i ) . 

Therefore, the probability of mission completion for the 
entire MAGTF is the product of the probabilities that the 
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individual secondary reparables exceed a level, k±, given 
spares (or equivalently experience at most 
failures). So, the maximum number of operational pacing 
items of a specific type at any given moment will be equal 
to the lowest level of its component secondary reparables. 
It is this minimum level that must exceed the specified 
level k since if the minimum level exceeds kj_, all levels 



will exceed kj_. Now, by the assumption that the secondary 



reparable failures are independent combined . with the fact 
that if one type of pacing item cannot complete the mission 
the entire MAGTF will not be able to complete the mission, 
the probability that the minimum level exceeds the specified 
level (i.e., p^(x^)) can be expressed as. 



which accounts for the functional form of the objective 
function. 

Similar to the first objective function, the second 
objective function also takes advantage of the assumptions 
of constant failure rates for each secondary reparable and 
independence among the failures of different secondary 
reparables; however, it incorporates a more realistic 
approach as to the nature of the failure rates. In the 
previous objective function, the failure rate A of each 
secondary reparable was assumed to be known (actually an 



L. 



.V 
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estimate computed using the maximum likelihood estimator of 



the failure rate must realistically be used). For this 
second objective function, a Bayesian approach is taken in 
which the exact value of ^ is explicitly recognized to be 
unknown. The observed data (i.e., times to failure) are 
used to update the prior information concerning each failure 
rate to form a posterior distribution as discussed in the 
chapter on data analysis. This posterior distribution is 
necessary to derive the predictive distribution and the 
functional form of Pj_(Xj_), since we are interested in 
computing the probabilty distribution of the number of 
failures of secondary reparable i which occur during some 
mission of length tj_. By analogy to the previous derivation 
of Pj_(x.j_) f it is desired to consider this number of 
failures because we want this quantity to be less than some 
specified level (i.e., a^+x^-k^) so as to insure the entire 
MAGTF completes the assigned mission. Assuming that the 
failure rate / ^ is given, the total number of failures of 
secondary reparables of type i in a mission of length 
will have a Poisson distribution with approximate mean 



But / i is not precisely known and has the following 
posterior density (as derived in the data analysis chapter) 




/’(Exactly f- l failures in (0,rj | / •) = — exp 

Jr 
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given the observed times to failure and prior information: 



A A 

A \n,- + a 



/(Ay I data) = 



ft* 4 - BY* 1 A .A 

= ..ini. lA X n i + «/ - 1 ovn + ft) 



I>y + a f ) 



exp 



where. 



T t total exposure time for secondary reparable i 

n i number of observations for secondary reparable 
i 

ay empirical Bayes estimate of shape parameter for Gamma 
( Ay ; ay , fii ) prior distribution 

A 

P ( empirical Bayes estimate of scale parameter for Gamma 
( Ay ; a. i , Pi ) prior distribution. 

This can be rewritten as, 

+ ' (T , + /?;) exp- ; ' <T ' + W 



/(Ay I data) = 



I . A y 

f (/Zy + CCy) 



Now, to calculate the predictive distribution, 

/’(Exactly jy failures in (0,rj) = f°° ^(Exactly failures in (0,fy]|Ay) /(Ay I data )d). i 

Jo 

which is upon substitution, 

f°°( k±£iL. cxp-^-VA ( < T f + h n ‘ t* 1 . ;»/ + Si - 1 exp _; -<< T ' + Wui, . 

Jo v r(», + £,) J 



Evaluation of the integral yields: 



/(Exactly^ failures in (0,ry]) = 



r(«y + ay +y~) / Ty *f /?y x + a ; 



( 



+ a i ^ 



a Ji 



Ji 



) • 



Jr + ay) rty/y -(- Ty + /?y 7 <?/y + Ty + /?y 

This is a negative binomial distribution with parameters 

T / + h 



P = 



fly/y + Ty + /?y 



and 
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r = Wy + ay . 
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Even though p^(x^) is defined as the probability that 
the stockage of x ^ spares for secondary reparable i allows 
the pacing item of which it is a component to exceed some 
minimum level (k^) , it has been shown that this is 
equivalent to the probability that the number of failures of 
secondary reparable i during mission t ^ is less than or 
equal to a^+x^-k^ (i.e., P( number of failures of secondary 

reparable i in a j_+ x j_-k ^) ) . Therefore, the second 

functional form of p^(x^) is: 

Q; ■+■ X: — k: A A a 

r(/7,- 4- Kj 4- tn) ( Tj 4- Pj ^ n i + *i ( a i l i 

°i( x i) = 2^ — ■ — — v r ) v 

m=0 mir ( n i + “/•) aft 4- t i 4- Pi aft + t ,• 4- />,• 

where, 

number of pacing items containing secondary reparable 
i 

minimum number of pacing items containing secondary 
reparable i required to be operational throughout 
given mission ( k^_< a^ ) 

number of observations ( data points ) of times between 
failure for secondary reparable i 

total observed exposure time for secondary reparable i 

mission duration during which secondary reparable i 
must operate 

empirical Bayes estimate of shape parameter for Gamma 
( A; ; a,- , Pi ) prior distribution 

empirical Bayes estimate of scale parameter for Gamma 
( /,• ; a,- , Pi ) prior distribution . 

The reasoning behind the functional form of equation (3.2) 

developed for the first p^(x^) still holds true for the 



a i 



k ; 



n, 



A 

a,- 



01 




37 



REPRODUCED AT GO VERNMFN I 



newly derived pj_(x^); therefore, the functional form of this 
objective function is still the product of the individual 
P±(x±) 's. 

The linear weight capacity restricts the total shipping 
weight (i.e., secondary reparable plus container) for all 
secondary reparables stocked to be less than or equal to the 
weight capacity of the ship or airplane. Similarly , the 
linear constraint on volume or space capacity requires that 
all secondary reparables stocked not exceed the volume 
available. The final constraint in PI indicates that there 
is an upper bound on the number of spares for each secondary 
reparable that can be stocked. This is a very real 
constraint because there are only a finite number of each 
type of secondary reparable at the wholesale inventory level 
from which the retail inventories used to support all 
exercises, deployments, and contingencies can be 
established. Thus, a combat service support planner must 
consider the requirements of the force as a whole so as not 
to commit too many secondary reparables to any one 
exercise . 

The number of constraints is dependent on the scenario. 
Any scenario that involves the airlift of the ground combat 
element and its supporting forces into the exercise or 
operation area, such as the short-term exercises described 
in the introduction, will have only one constraint. This is 
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the constraint on the allowable cargo load (ACL). Allowable 
cargo load is: 

The initial limiting factor in determining airlift 
requirements. . . . This is a weight limiting factor 

expressed in terms of STON [short tons] which cannot be 
exceeded in planning aircraft loads. The ACL is 
established by the Air Force for each type aircraft for a 
specific mission. [Ref. 22:pp. 4-9] 

Since the formulation of problem PI now reduces to just one 

constraint, it is recognizable as a nonlinear knapsack 

problem due to the nonlinear objective function. The 

individual secondary reparables represent the items or 

commodities that must fit into a knapsack while the ACL 

constraint corresponds to the knapsack's weight or volume 

capacity. 

The other two scenarios described in the introduction 
both involve two constraints. These two constraints 
represent the weight and volume capacities of the shipping 
assigned to perform the sealift. In the case of the MAUs 
deploying to the Mediterranean Sea, the secondary reparables 
are normally stored aboard one of the ships of the 
amphibious task force such as a LHA, LPH, LKA, LSD, or LST . 
Thus, the constraints reflect the weight and volume capacity 
of the ship's cargo hold or compartment which has been 
designated for secondary reparable storage. On the other 
hand, in considering secondary reparable support of the MPF 
MAB, the secondary reparables are spre ad- loaded among 
several of the MPF squadron ships; therefore, the weight and 
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volume constraints are actually the combined weight and 
volume capacities of the cargo holds of the specified MPF 
squadron ships. Since both of these sealift scenarios 
involve two contraints, problem PI can no longer be regarded 
as a nonlinear knapsack problem, but it is still an integer 
nonlinear programming problem. 

Regardless of the number of constraints or whether 
equation (3.1) or (3.3) is used as the functional form for 
pj_(Xj_), the objective function is a complex, nonlinear 
function of the decision variable. For this reason, problem 
PI will be transformed into a binary linear programming 
problem. The following section explains this transformation 
and the new formulation. 

B. TRANSFORMATION TO A BINARY LINEAR PROGRAM 

It can be seen from the final constraint of problem PI 
that the decision variable, x^, is constrained to a set of 
discrete values. Therefore the objective function will 
first be separated using the properties of logarithms and 
then the transformation technique for converting discrete 
variables to binary variables will be used to reformulate 
the problem [Ref. 23:p. 12]. Taking the natural logarithm 
of the objective function yields: 




Now, the final constraint in PI can be expressed as. 



X{ 6 S{ = {Sj-g, , ... , 5/y.} — {0,1, ... ,W/} 

and this is equivalent to the constraint set [Ref. 

12 ] : 

u i 

v 

x i= ^ifij 
i= 0 

£** = > 

/=0 



23 : p 



X;: e {0.1} , j- 0,1 U[ . 



Applying the separation technique and transformation to 
problem PI results in the new formulation: 



INDICES: 
/= 1, ... ,.Y 

j = 0 u { 

data: 

fij 

W V 

IV 



v iJ 



U: 



secondary reparable or component 
number of secondary reparables i 

ln(p ± ( j ) ) 

total shipping weight for j spares of type i 
[lbs] 

weight capacity of airplane or ship [lbs] 

total volume of j spare of typei[cuft] 

volume capacity of airplane or ship [cuft] 

maximum allowable number of spares for 
component i 



DECISION VARIABLE: 



x ij 



(.0 othe 



spares of type i are selected 
otherwise 
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FORMULATION: 

S u i 

max V V f ijX .. 

<=1 y=0 



(P 2) 



s.t. 



-V 

E Z-Vff s w 

(=1 7=0 



y y % s k 



i=l 7=0 



u i 

E^=‘ v/ 

y =0 

A-y e (0,1} 

Assuming known failure rates, the coefficients, f±j' i n the 
objective function are: 



fii = In 



c a i 

V . 

/_J 

m = 0 



(a/.;,) 



m 



A77: 



exp 



A 

-Wi 



( 3 . 4 ) 



where, 

<*i number of pacing items containing secondary reparable 
i 

k{ minimum number of pacing items containing secondary 
reparable i required to be operational throughout 
given mission (k^<a^) 

A 

M.L.E. of failure rate for secondary reparable i 

l i mission duration during which secondary reparable i 
must operate. 



If perfect knowledge about the failure rates is not assumed 
and the Bayesian approach is adopted the coefficients are: 



42 



/ 



r a i +/- k i 



4 = 1 n 



y 



n«/ + a,- + m) 



(-^ )"*'■( 

a.Y: + T, -f /?; 



A 



*- m 



=0 w! ffo-rog 



a i r i 



a i l i 



T . 

"L 



Pi 



m 



( 3 . 5 ) 



where , 

a,- number of pacing items containing secondary reparable 
i 

k i minimum number of pacing items containing secondary 
reparable i required to be operational throughout 
given mission ( k^£a^ ) 

n i number of observations (data points) of times between 
failure for secondary reparable i 

T i total observed exposure time for secondary reparable i 

i[ mission duration during which secondary reparable i 
must operate 

a ; - empirical Bayes estimate of shape parameter for Gamma 
( X ±f & ±, p ±) prior distribution 

A 

/?• empirical Bayes estimate of scale parameter for Gamma 
( j_; 2 j_ , /?j_) prior distribution. 

It is of interest to note that for situations in which 
the airlift model is applicable, problem P2 is just a 
special case of the generalized assignment problem. In the 
formulation of the generalized assignment problem 
[Ref. 24 : pp . 345-346] , each agent (or man) can perform or be 
assigned more than one task (or job) provided that the 
resource available to the agent is not exceeded; whereas, in 
the classical assignment problem [Ref. 23:pp. 61-62], each 

agent can be assigned to only one task. In both 
formulations, each task must be assigned to exactly one of 
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the agents. The airlift model can be interpreted as a 
special case of this generalized assignment problem in which 
one agent must perform all N tasks, however, he has (u^ + 1) 
alternative ways to accomplish each task. 

Problem P2 is equivalent to PI and even though the 
number of variables and constraints have increased, it is 
the preferred formulation because it eliminates the complex 
nonlinear objective function of PI. The coefficients fj_j 
can be directly computed using recursive formulas for the 
Poisson or negative binomial cumulative distribution 
functions; whereas, Pj_(Xj_) could not since it was itself a 
function of the decision variable. Additionally, because of 
the assumption of the independence of failures, other 
failure distributions can be readily introduced into this 
linear formulation. Furthermore, the formulation is general 
enough so that additional constraints can be incorporated as 
required and as will be seen in the following chapter, it 
lends itself easily to programming in an algebraic modeling 
language. 
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IV. SOLUTION METHODOLOGY AND IMPLEMENTATION 



The purpose of this study is to determine an optimal 
stockage level of spare secondary reparables for each of the 
three different scenarios described in the introduction. 
Of course "optimality" is with reference to the failure 
model adopted and the relevance of the data used to estimate 
its parameters. The determination of the stockage levels 
is the subject of this chapter. The binary linear nature of 
the formulation presented in the preceding chapter suggests 
that several different methods of solution are available to 
determine the optimal stockage policy for a given scenario. 
Each alternative method will be discussed and justification 
for the method chosen is provided in addition to its 
implementation . 

A. SOLUTION APPROACHES AND METHODOLOGY 

The mode of deployment in a particular scenario 
determines the number of constraints. The initial solution 
approach to the scenario in which the airlift model was 
applicable utilized the technique of Lagrangian relaxation. 
Although it is standard procedure to establish bounds for 
integer programs (such as problem P2) using a linear 
programming ( LP ) relaxation, it was initially desired to 
avoid dependency on commercial LP solvers not locally 
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available to the combat service support planners of the 
Second Force Service Support Group at Camp Lejeune, North 
Carolina. Therefore, Lagrangian relaxation was adopted 
because it lends itself nicely to programming in readily 
available FORTRAN 77. It must be emphasized that the 
structure of problem P2 did not demand that this specific 
r solution technique be applied. In fact, the LP relaxation 

of problem P2 could have been solved using an LP solver. 
Thus, Lagrangian relaxation was utilized not out of 
necessity, but primarily because it offered a realistic 
alternative to the use of a commercial LP optimizer. 

Briefly stated, the technique of Lagrangian relaxation, 
as outlined by Fisher [Ref. 25:pp. 10-21], involves moving 

complicating constraints into the objective function using 
the product of the Lagrangian multiplier and the constraint 
violation as a penalty term. In particular, examination of 
problem P2 in the context of the airlift scenario results in 
the following binary linear formulation (actually a 
generalized assignment formulation): 

INDICES: 

/_ 1 ... n secondary reparables or components 
y = 0, number of secondary reparables i 
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DATA : 



Jij 


ln(p ± ( j ) ) 






M V 


total shipping weight for j 
[lbs} 


spares of type i 


ACL 


allowable cargo load of airplane 


[lbs] 


u i 


maximum allowable number 
component i 


of 


spares for 


DECISION 


VARIABLE: 







x .. _ fl if j spares of type i are selected 
lJ [0 otherwise 



FORMULATION: 

-V u i 

v YjYj fir x ij 
i= 1 7=0 



max 



.V u i 

>./■. / / < ACL 

/= 1 7=0 



(P3) 




V/ 



7=0 

Xije { 0 , 1 } 

The only complicating constraint is the allowable cargo load 
(ACL) constraint. The relaxed formulation with s_cal ar 
Lagrangian multiplier, A , consists of N separable, muJJzijrte 
choice problems which are easy to solve. On the other hand, 
for those situations in which both the weight and. volume 
constraints of problem P2 must be considered (i.e., 
Mediterranean MAU or MPF MAB scenarios), the Lagrangian 
multiplier, A , is now vector-valued (i.e., A =( A]_, A 2) since 
two constraints must be relaxed, and this is no longer a 
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trivial problem to solve. The primary reason for this 
complication is that the solution approach now involves a 
two-dimensional search which is more difficult to program 
than the simple bisection search associated with the single, 
scalar Lagrangian multiplier. The problem is not impossible 
to solve and methods exist, such as a coordinate search 
method suggested by DeWolfe, Stevens, and Wood [Ref. 26: 
pp.8-10], which appears to offer good results for this type 
of problem. 

As additional constraints (such as budget, holding 
(storage) cost, or shortage (stockout) cost constraints) are 
added, the programming difficulty associated with the search 
procedure of the Lagrangian function becomes prohibitive. 
Therefore, to accommodate the possible inclusion of 
additional constraints in the formulation of problem P2, it 
was decided to utilize the Generalized Algebraic Modeling 
System (GAMS) in conjunction with the Zero/One Optimization 
Methods (ZOOM) mixed integer program solver. The use of 
this optimization package not only avoids the restrictive 
programming considerations discussed above, but it can also 
incorporate the PEB as an initial step, prior to solution of 
P2 or P3 . 

The approach to the PEB is to maximize the marginal 
likelihood function associated with each secondary reparable 
i with respect to and This marginal likelihood 
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function is a nonlinear function of the two decision 



variables. The formulation for this maximization problem 
is : 

INDICES : 

,N secondary reparable or component 
m = 1 Mj source of data (e.g., tank serial number 508282) 



DATA: 

number of observations of component i failures 
for data source m 

total observed exposure time for component i of 
data source m. 

upper bound on the shape parameter of the Gamma 
prior distribution for component i. 

Pi upper bound on the scale parameter of the Gamma 

prior distribution for component i. 



n, 



im 



im 



DECISION VARIABLES: 

A 

a i shape parameter for Gamma prior distribution 

associated with compenent i. 

A 

Pi scale parameter for Gamma prior distribution 

associated with component i. 



FORMULATION: 

Ml 

max 



r(n im t ocy) 

'W n«,) 




* tm 



l im Pi 




(M) 



s.t. 0 <a i <a. i 



0 < Pi < Pi 



Standard maximization by taking partial derivatives with 
respect to and does not yield a satisfactory closed 
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form solution which can be practically implemented. Thus, 
this nonlinear maximization problem is also solved using 
GAMS in conjunction with the nonlinear program solver, 
MINOS [Ref. 27] . These values are then incorporated into 
the optimization of problem P2. In summary, under the 
current solution approach, the entire model can be coded for 
computer processing on either a mainframe system or personal 
computer and solved using a single modeling system (i.e., 
GAMS) and its compatible solvers. This approach is 
sufficiently flexible to incorporate additional constraints 
or decision variables. 

B. IMPLEMENTATION 

Since the decision was made to determine the optimal 
stockage policy for all scenarios using GAMS, no specifics 
of the original solution approach involving Lagrangian 
relaxation will be mentioned other than the fact that it is 
implemented by means of a FORTRAN 77 computer program. The 
FORTRAN program is based upon a similar one developed by 
DeWolfe [Ref. 17:pp. 23-24]. A copy of the source code 
developed for this study is provided in Appendix A and 
examples of its two required input files and of the single 
output file are given in Appendices B and C, respectively. 

The solution approach adopted centers around the coding 
of problems P2 and P4 for computer processing using the GAMS 
language. The reader is referred to the excellent tutorial 
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by Rosenthal [Ref. 28] for the specifics of this versatile 
mathematical modeling language. Upon coding, the problems 
are solved using the GAMS compatible solvers of MINOS and 
ZOOM for nonlinear programs and mixed integer programs 
respectively. Specifically, the source code provided in 
Appendix D incorporates the successive solve feature of GAMS 
to first solve the nonlinear maximization problem 
represented by problem P4. The results of this maximization 

A A 

provide the estimates of the shape ( a^) and scale ( p 
parameters of the prior distributions which are required 
before the second binary linear maximization (P2) can be 
performed. The final result is the optimal stockage level 
of spare secondary reparables for the given scenario. Even 
though the source code in Appendix D is specific to the 
airlift model ( P3 ) , GAMS makes it very easy to incorporate 
additional constraints. For example, if a volume constraint 
is also applicable, a single GAMS PARAMETER statement to 
accomodate the volume data and EQUATION statement to reflect 
the volume constraint are the only additions necessary. 

The source code provided is compatible with either 
mainframe or personal computer usage. When run on the IBM 
3033AP mainframe computer, GAMS version 2.05 was utilized in 
conjunction with MINOS Version 5.1 and ZOOM version 2.1 and 
on the Zenith Z-248 personal computer (equipped with 80287-8 
math coprocessor and 640 Kbyte memory expansion board), PC 
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GAMS Version 2.04 was utilized in conjunction with MINOS 
Version 5.0 and ZOOM Version 2.1. 



c 
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V. RESULTS 



This chapter summarizes results obtained using actual 
operational maintenance data to determine an optimal level 
of spare secondary reparables for a Marine Amphibious Unit 
(MAU) in a realistic exercise scenario. Verification of the 
model is accomplished using simulated failure data. 
Estimation of the parameters of the prior distribution by 
an analysis of the data is discussed and demonstrated; this 
realizes computationally the PEB procedure. Finally, 
computational experience for both the mainframe and Zenith- 
248 personal computer is briefly discussed. 

A. MODEL TESTING 

The data used in the development and testing of all 
models was provided by the Operations Section of the 
Supported Activities Supply System (SASSY) Management Unit 
in coordination with the Maintenance Information Systems 
Management Office of the Second Force Service Support Group 
at Camp LeJeune, North Carolina. This data was extracted 
from the ERO History File of the Marine Corps Integrated 
Maintenance Management System and analyzed according to the 
procedure presented in Chapter II. This analysis provided 
the basic parameters for equation (3.5) (or (3.4)) from 
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which the objective function coefficients were generated 
using GAMS. 

The airlift model ( P3 ) was used to determine the optimal 
stockage policy for twenty-five different secondary 
reparables associated with the pacing items of a Marine 
Amphibious Unit participating in a realistic short-term 
exercise. This problem generates 190 binary variables and 
26 constraints. The recommended stockage level for the 
secondary reparables associated with the three pacing items 
(i.e., M60A1 Tank, LTP7A1 Amtrack, and M101A1 105mm 
Howitzer) in this exercise is presented in Table 2. Table 2 
displays the levels computed using both the GAMS/ MINOS/ZOOM 
package and the Lagrangian relaxation technique. 

For this exercise scenario 99.89% of the allowable cargo 
load was utilized and the objective function value was 
within 0.07% of the upper bound established during the 
branch and bound procedure used by the ZOOM solver. The 
probability of mission completion given the stockage level 
was 0.4149 and this highlights one of the advantages 
associated with using this model. This low probability 
implies that the allowable cargo load must be increased. The 
only way to accomplish this is to increase the number of 
cargo planes assigned to this exercise. This provides 
Marine Corps planners with a tangible basis on which to 
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TABLE 2 



MAU AIRLIFT SCENARIO SPARE STOCKAGE LEVELS 

SPARE STOCKAGE LEVEL 

GAMS/ LAGRANGIAN 
NSN END ITEM ZOOM RELAXATION 



2815001245387 


M60A1 


2 


2 


2520002241867 


M60A1 


3 


3 


2815002395819 


M60A1 


2 


2 


2520001549632 


M60A1 


2 


2 


2530000886657 


M60A1 


2 


3 


2520001184942 


M60A1 


2 


2 


2920010135802 


M60A1 


1 


1 


1010010703803 


M60A1 


1 


1 


1025019388105 


M60A1 


4 


5 


1010010720397 


M60A1 


1 


1 


1025010708993 


M60A1 


1 


1 


2815004303480 


LVTP7A1 


3 


4 


2520003973384 


LVTP7A1 


5 


5 


2520001443385 


LVTP7A1 


4 


4 


5805014593410 


LVTP7A1 


4 


5 


2520008949532 


LVTP7A1 


2 


2 


2910011714636 


LVTP7A1 


2 


3 


2530011509757 


LVTP7A1 


2 


2 


5865012207848 


LVTP7A1 


3 


3 


4320012035660 


LVTP7A1 


2 


2 


2920002317276 


LVTP7A1 


1 


1 


2530000886650 


LVTP7A1 


2 


2 


1005011855059 


M101A1 


5 


4 


1005011086434 


M101A1 


5 


5 


1005011457709 


M101A1 


5 


5 
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justify their requests for additional airlift assets from 
the Military Airlift Command of the U.S. Air Force. Table 3 
summarizes the results for both solution approaches. 

TABLE 3 

SUMMARY OF RESULTS FOR MAU AIRLIFT SCENARIO 

Objective Probability Deviation 

Solution Function of Mission Percentage from 

r Technique Value Completion of ACL Used Optimality 

GAMS /MINOS/ 

ZOOM -0.8797 0.4149 99.89% 0.07% 

Lagrangian 

Relaxation -0.8820 0.4139 100% 0.47% 

( w/Heuristic ) 

B. MODEL VERIFICATION 

Verification of the modeling approach was conducted by 
using the model to verify the fact that the negative 
binomial predictive function reduces to the Poisson-Gamma 
function as the number of failures observed (and accordingly 
the total exposure time) approaches infinity. Thus, the 
Poisson function provides the limiting case or best possible 
stockage level attainable assuming that the failure rates 
for each secondary reparable are precisely known. 

To accomplish this verification, the MAU airlift model 
was run using the Poisson objective function (3.4) with 
fixed failure rates. The resulting stockage levels 
represent the best stockage policy that can be achieved. 
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Then, using those failure rates, twenty-five different sets 
of simulated failure data were generated for each of 10, 25, 
50, 100, 200 and 500 failures. These simulated failure 
times were randomly generated using the probability integral 
transform to obtain exponential failure times from uniform 
(0,1) deviates provided by the LLRANDOMII Random Number 
Generation Package [Ref. 29] . The Bayesian analysis of this 
data was performed resulting in the negative binomial 
predictive distribution and the objective function, equation 
(3.5). A noninformative Gamma prior distribution with a =1 
and (j =1 was employed in this process. Stockage levels for 
this method are displayed in Table 4. Note that about 100 
failure observations are needed before the stockage levels 
approach those computed using the actual failure rates 
(i.e., the optimal stockage policy given perfect 
information). Figure 1 is a plot of the interquartile range 
(IQR) of the probabiltiy of mission completion associated 
with each failure sample size. This graph demonstrates that 
as the number of observed failures increases, the 
variability (represented by the IQR) in the probability of 
mission completion decreases. 

In addition to the above verification, the fact that the 
model and solution approach responded predictably to a wide 
range of different failure observations and exposure times 
indicates that the model and approach tend to be robust. 
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TABLE 4 



MAU AIRLIFT SCENARIO SPARE STOCKAGE LEVELS 
SPARE STOCKAGE LEVEL 



NSN 



NEGATIVE BINOMIAL 

10 25 50 100 200 500 POISSON 



2815001245387 

2520002241867 

2815002395819 

2520001549632 

2530000886657 

2520001184942 

2920010135802 

1010010703803 

1025019388105 

1010010720397 

1025010708993 

2815004303480 

2520003973384 

2520001443385 

5805014593410 

2520008949532 

2910011714636 

2530011509757 

5865012207848 

4320012035660 

2920002317276 

2530000886650 

1005011855059 

1005011086434 

1005011457709 



3 

4 
2 
2 
1 
1 
1 
0 

4 
1 
1 

5 
5 
3 
5 
2 
1 
2 
3 

3 
1 
1 

4 

5 
5 



3 

4 
2 
2 
1 
1 
1 
0 
4 
1 
0 

4 

5 

4 

5 
2 
2 
2 
2 
2 
1 
1 

4 

5 
5 



3 

4 
2 
2 
1 
1 
1 
1 
4 
1 
0 

4 

5 

4 

5 
2 
2 
2 
2 
2 
1 
1 
3 
5 
5 



2 

3 

2 

3 
2 
2 
1 
1 

4 
1 
1 

4 

5 

4 

5 
2 
2 
2 

3 
2 
1 
1 

4 

5 
5 



2 

3 

2 

3 
2 
2 
1 
1 

4 
1 
1 

4 

5 

4 

5 
2 
2 
2 

3 
2 
1 
1 

4 

5 
5 



2 

3 

2 

3 
2 
1 
1 
1 

4 
1 
1 

4 

5 

4 

5 
2 
2 
2 

3 
2 
1 
2 

4 

5 
5 



58 



/ 






PROBABILITY OF MISSION COMPLETION 
0.35 0.40 0.45 0.50 0.55 0.60 




Figure 1. Variability in P( Mission Completion) with 
Number of Observed Failures 



59 



RETROD U CCD AT GOVERN MFD I 



C. BAYESIAN AND EMPIRICAL BAYES ANALYSIS 

Bayesian analysis and particularly PEB suggest 
themselves when the failure rates of different units (e.g., 
tank engines on different tanks or amtrack transmissions on 
different amtracks) are sufficiently alike, but individual 
experience is small, so that "borrowing strength" by a 
r suitable pooling process can improve individual estimates. 

In this study, the failure rate associated with each 
secondary reparable ( /. ^ ) was assumed to follow a Gamma 

prior distribution. This prior distribution represents the 
probability that the unknown value of the failure rate lies 
in various regions of the parameter space [Ref. 20 :p. 260] . 
In other words, incorporating the prior distribution enables 
a tighter bound on the unknown value of the failure rate to 
be obtained than would be possible by using the data 
pertaining to that rate alone. The primary disadvantage of 
the ordinary or classical Bayesian approach is that the 
prior distribution's parameters (i.e., oc^ and are 

assumed to be known. This deficiency is overcome in this 
study by employing a PEB approach, in which the entire data 
set ( {n im> t im ). m = 1 , ... , Mi ) associated with each secondary 
reparable i is used to estimate the parameters of the prior 
distribution. These estimates are then used in the standard 
Bayesian analysis to compute the posterior and predictive 
distributions. The results are unavoidably influenced by 
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the particular choice of prior used (here a Gamma). 
Sentitivity tests have not been conducted, but indication is 
that the procedure is rather insensitive. 

A specific indication of this insensitivity is seen by 
comparing the airlift model stockage levels obtained when 
minimal prior information is included to the levels obtained 
when the empirical Bayes estimates of the prior distribution 
are included. The results of this comparison are summarized 
in Table 5. 

The levels in the column labeled "NONINFORMATIVE PRIOR" 
were obtained by running the airlift model with the randomly 
generated failure times (described in the previous section) 
associated with the number of failures sample size of ten. 
A "diffuse" or noninf ormat ive prior distribution 

(i.e., y. i = l ) was used. A diffuse prior distribution is 
"informationless in a relative sense and means only that it 
is diffuse relative to the sample information" . [Ref . 30:p. 

198] The same airlift model was then run using empirical 
Bayes estimates for the prior parameters and these levels 
appear in the column labeled "INFORMATIVE PRIOR". These 
levels are only slightly closer to the Poisson levels than 
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TABLE 5 



EMPIRICAL BAYES ESTIMATION 



STOCKAGE LEVELS 

NEGATIVE BINOMIAL POISSON 

NONINFORMATIVE INFORMATIVE 
NSN PRIOR PRIOR 



2815001245387 3 
2520002241867 4 
2815002395819 2 
2520001549632 2 
2530000886657 1 
2520001184942 1 
2920010135802 1 
1010010703803 0 
1025019388105 4 
1010010720397 1 
1025010708993 1 
2815004303480 5 
2520003973384 5 
2520001443385 3 
5805014593410 5 
2520008949532 2 
2910011714636 1 
2530011509757 2 
5865012207848 3 
4320012035660 3 
2920002317276 1 
2530000886650 1 
1005011855059 4 
1005011086434 5 
1005011457709 5 



2 

4 

2 

2 

1 

1 

1 

0 

4 
1 
1 

5 
5 

4 

5 
2 
1 
2 
3 

3 
1 
1 

4 

5 
5 



2 

3 

2 

3 
2 
1 
1 
1 

4 
1 
1 

4 

5 

4 

5 
2 
2 
2 

3 
2 
1 
2 

4 

5 
5 
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are those computed when minimal prior information is 
considered. 

Even though these stockage levels are only slightly 
closer to the Poisson levels for this specific instance, by 
simply using the available data for 10 failures to compute 
empirical Bayes estimates of the prior distribution 
parameters, spare stockage levels comparable to those 
obtained in the previous section using failure sample sizes 
of between 50 to 100 can be achieved. This is particularly 
significant since we have no control over the number of 
failures observed for each secondary reparable or their 
exposure times and would not usually observe this number of 
failures (i.e., 50-100) during the reporting period covered 
by the ERO History File. Also, since operational 
reliability data may not be available for the exact location 
where the MAGTF is deploying, pooling the data from sources 
which have operated in similar environments has the 
advantage of providing more realistic failures rates than if 
this prior information were not considered. 

D. COMPUTATIONAL EXPERIENCE 

So as to keep the implementation and solution of the 
models developed in this study as general and flexible as 
possible existing software packages (i.e., STATGRAPHICS , 
GAMS, etc.) were utilized. The goal was to obtain solutions 
in a relatively user-friendly personal computer environment 
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requiring a minimum of operator interaction. Therefore, 
computational efficiency in terms of computer run time (or 
algorithmic complexity) was not of primary concern and is 
included here solely for completeness. No attempt was made 
to employ special algorithms that exploit problem structure. 
This is not to downplay the importance of this factor, but 
rather to emphasize that during this study the primary 
concern was to obtain solutions to real-world problems as 
easily and efficiently as possible with a minimum of 
programming effort. 

Running time for the airlift model on the IBM 3033AP 
using the GAMS/MINOS/ZOOM package was 3.9 seconds in 
contrast to the 176 seconds obtained when run on the Zenith 
Z-248 PC. It is possible on the IBM 3033AP to combine the 
PEB calculations with the maximization of problem P3 by 
using the successive solve feature of GAMS; however, this 
cannot be done on the Zenith Z-248 PC since internal 
parameters in the PC version of GAMS are exceeded. 
Therefore, each of the empirical Bayes maximizations must be 
performed separately when using this PC. The resulting 
estimates have to be manually input into the appropriate 
GAMS PARAMETER statements before maximization of problem P3 
can be accomplished. 

Even though the Lagrangian relaxation procedure with 
heuristic was not adopted as the solution technique for the 
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reasons discussed in Chapter IV, it provides much faster 
solution times to the airlift model ( P3 ) on both the IBM 
3033AP and Zenith Z-248 PC. The Lagrangian procedure, coded 
in FORTRAN 77, provides a running time of 0.32 seconds on 
the IBM 3033AP and about 30 seconds on the Zenith Z-248 PC 
where it was coded using Microsoft FORTRAN (Version 3.2). 

Finally, comparison of the stockage levels obtained 
using GAMS/MINOS/ZOOM with those obtained using LP 
relaxation for each of the 150 simulated failure time cases 
revealed that in every case the LP relaxation provided 
identical stockage levels with the exception that it 
fractionated for exactly one of the twenty-five secondary 
reparables. So, if a general policy was adopted to round 
down to the lower of the fractionated stockage levels ( so as 
not to violate the ACL contrast), the LP relaxation would 
provide comparable results to those obtained using the mixed 
integer program solver, ZOOM. This is an important result 
because the run times associated with the LP relaxation were 
significantly less than those using ZOOM. In particular, 
if the rounding down policy is applied to the LP relaxation 
of the airlift model in this study, the resulting integer 
stockage level provides an objective function value which is 
within 1.87% of the LP upper bound or 1.73% of the upper 
bound established by the branch and bound procedure utilized 
by ZOOM. This is not quite as good as the result obtained 
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directly from ZOOM (within 0.07% of optimality), however, 
for those situations in which computer run time may be a 
consideration (e.g., as problem size increases) or in which 
very tight bounds may not be required, a very good solution 
can be obtained from the much simpler and less costly LP 
relaxation. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



In this thesis, the problem of determining spare 
stockage levels for deployable MAGTFs was 
solved by 

(a) stochastically modeling demand utilizing operational 
data from the Marine Corps Integrated Maintenance 
Management System, and then 

(b) using mathematical programming to optimize a 
suitable, operationally relevant objective function. 

To accomplish this, a new measure of effectiveness was 

defined which provides the operational commander on the 

ground with an objective indicator of the overall impact of 

the stockage policy on the probability of mission 

completion. This measure of effectiveness was evaluated as 

the end result of a probabilistic modeling step, a 

statistical data analysis, and a subsequent optimization. 

The generality of the model formulation makes it 
adaptable to any of the three scenarios of present interest 
to the Marine Corps. The use of a PC-compatible algebraic 
modeling language (GAMS) makes it particularly easy to add 
or delete constraints as required. Specifically, the model 
and solution approach were applied to a realistic exercise 
scenario and an essentially optimal stockage policy was 
determined which maximized the probability of mission 
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completion within the constraints imposed by the deployment 
mode. 

During the analysis of the MIMMS data base (i.e.,ERO 
History File) numerous problems and inconsistencies were 
encountered. As a result, the estimates of the failure 
rates for some of the secondary reparables are higher than 
would be the case if all the potential times to failure 
could be extracted from the data base. To overcome this 
discrepancy, it is highly recommended that the quarterly 
update procedure described in Chapter II be conducted at the 
Force Service Support Group level. If this cannot be done, 
then at the very least the meter readings associated with 
all pacing items should be recorded and input into MIMMS 
prior to each quarterly update. This would insure that 
beginning and ending meter readings are available for each 
serialized pacing item during the reporting period covered 
by the ERO History File and this would prevent the loss of 
valuable exposure time data. 

The limited strategic lift assets available to deploying 
Marine Corps units demands that effective utilization be 
made of these resources. Since only a small percentage of 
these assets can be allocated for spare parts inventories, 
emphasis must be placed on stocking those parts which would 
have the greatest impact on the completion of the assigned 
mission. The reliability-based model developed in this 
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study accomplishes this objective by maximizing the 
probability of mission completion within the constraints of 
the mode of deployment. 
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APPENDIX A 



LISTING OF SOURCE CODE FOR LANGRANGIAN RELAXATION 



PROGRAM RELAX2 

CHARACTER NSNX13,N0MENX7 

INTEGER NCOMP, NEQUIP , NSPARE, MAXSPR, IDNUM, BSP ARE, SPBND, NDENS, 
XN0PER,NXP0SE 

REAL X8 P,W, MAXWT, DINF, OPTDEV, ZL0W,ZUP, PSUC, HI LAMD,WT,MTBF, 

XDURMIS,ZHLOW,XPOSE, ALPHA, BETA 

PARAMETER ( NC0MP=25 , NEQUIP=3, MAXSPR=1 0 , MAXWT=4440 0 . 0D0 , 

XDINF=1 . 0D20 ) 

DIMENSION NSN (NCOMP), NOMEN (NCOMP),NSPARE( NCOMP ), IDNUM ( NCOMP ) , 
XP(NCOM?,0 -.MAXSPR) ,W(NCOMP, 0 :MAXSPR) , BSPARE( NCOMP) , SPBND( NCOMP) , 
XNDENS ( NCOMP), N OP ER( NCOMP), XP0SE( NCOMP ), NXP0SE( NCOMP) , 

XDURMIS( NCOMP), ALPHA (NCOMP), BETA (NCOMP) 

CALL I NI T2 ( DINF, MAXSPR, NEQUIP , MAXWT, P , W, NSN, NOMEN, HI L AMD » NSPARE, 
XSPBND, NCOMP, IDNUM, NDENS, NOP ER, XPOSE, NXPOSE, DURM IS, ALPHA, BETA) 

CALL BOUND(DINF,P,W, NCOMP, MAXSPR, SPBND, MAXWT, HI LAMD, NSPARE, ZLOW, 
XZUP,WT,BSPARE) 

CALL HERISK DINF, P,W, NCOMP, MAXSPR, ZHLOW,WT, MAXWT, SPBND, BSPARE, 

XL ARRAY ) 

WRITE (10,50) 



50 FORMAT ( ' 1 '/ ' 0 ' , » 

x' » ) 

OPTDEV = -1.0D0 X 100. 0D0 x (1.0D0 - (ZHLOW/ZU P)) 



WRITE (10,100) ZHLOW 

100 FORMAT ( IX, 'OBJECTIVE FUNCTION VALUE IS :',D17.10) 

WRITE (10,150) OPTDEV 

150 FORMAT ('O’. 'THIS SOLUTION IS WITHIN ',D17.10,' % OF THE OPTIMAL') 
PRINT X, »ZUP=’,ZUP, 'ZLOW*', ZLOW, ’ ZHLOW= ', ZHLOW 
PSUC = EX P( ZHLOW ) 

WRITE (10,175) PSUC 

175 FORMAT ( ' 0 ' , ' P( MISSI ON SUCCESS) = ',F7.5) 

WRITE (10,150) 

ISO FORMAT ( '0' , ' ' , 

x' ' ) 

WRITE (10,200) 

200 FORMAT ( ’ 0 ’ , 5X, ’ NSN ’ , 10X, ’ END ITEM ', 6X, ' SPARE LEVEL’) 

DO 250 1=1, NCOMP 

WRITE (10,300) NSN( I ) , NOMEN( IDNUM( I ) ) , BSPARE( I ) 

250 CONTINUE 

300 FORMAT ( ' 0 ’ , A13 , 7X, A7 , 7X, 14) 

PRINT x, 'WEIGHT USED= ', WT, ' MAXWT= MAXWT 

STOP 

END 
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10 

20 

21 

22 

50 

25 



125 

100 

75 



SUBROUTINE INIT2C DINF, MAXSPR , NEQUIP , MAXWT , P . W, NSN , NOMEN , HI LAMD, 
XNSPARE,SPBND, NCOMP, I DNUM,NDENS,NOPER,XPOSE,NXPCSE, DURMIS, ALPHA, 
XBETA) 

INTEGER NCOMP 

REAL*8 MAXWT 

INTEGER IDNUM(NCOMP),NDENS(NCOMP) , NOPERC NCOMP ) , LMT, MAXSPR, 
XNSPAREC NCOMP), SPBNDC NCOMP) , NEQUIP, NXPOSE( NCOMP) 

REALMS WC NCOMP, 0 : MAXSPR ), DURMIS ( NCOMP ) , DINF, PI, Q, BET A (NCOMP), 
XPTEMP1,P(NCOMP,0: MAXSPR) , HI LAMD, HLTEMP, CUMSUM, XPOSEC NCOMP ),PP, 
XALPHACNCOMP) 

CHARACTER NSNX13(NC0MP),N0MENX7(NC0MP) 

READ (1,10) ( NSN( I ) , IDNUMC I ) , XPOSEC I) , NXPOSE( I ) , W( I , 1 ) , SPBNDC I ) , 
XALPHA(I),BETA(I), 1=1, NCOMP) 

FORMAT (A13,2X,I2,2X,F9.1,2X,I4,2X,F6.1,2X,I3,2X,F7.4,2X,F8.2) 

READ (2,20) (NOMEN(I),NDENS(I),NOPER(I), DURMIS ( I ) , 1=1, NEQUIP) 

FORMAT (A7,4X,I3,4X,I3,4X,F6 .1) 

DO 21 1=1, NCOMP 

PRINT x,NSN(I),IDNUM(I),XPOSE(I),NXPOSE(I),W(I,l),SPBNDCI), 
XALPHA(I),BETA(I) 

CONTINUE 

DO 22 1=1, NEQUIP 

PRINT x, NOMEN ( I ) , NDENSC I ) , NOPERC I ) , DURMIS ( I ) 

CONTINUE 



DO 25 1=1, NCOMP 
NSPAF.ECI) = 0 
DO 50 J=0, MAXSPR 
P ( I , J ) = DINF 

IF (J .NE. 1) W(I,J) = DINF 
CONTINUE 
CONTINUE 

DO 75 1=1, NCOMP 

PI = (XPOSE(I)+BETA(I))/(NDENS(IDNUM(I))x 
XDURMISC IDNUMC I) ) + XPOSEC I) + BETAC I ) ) 

Q = 1 . 0D0 - PI 

PP = PI XX C NXPOSEC I ) + ALPHA C I ) ) 

LMT = NDENSC IDNUMC I ) ) - NOPERC IDNUMC I) ) 

DO 100 J=0 , SPBNDC I ) 

IF CJ .NE. 1) W(I,J) = JXWC I , 1 ) 

CUMSUM = 0.0D0 
DO 125 K=0 , ( LMT + J) 

IF (K .EQ. 0) THEN 
PTEMP1 = 1 . 0 DO 
ELSE 

PTEMP1 = PTEMP1XCCCCNXPOSECI)+ALPHA(I))+(K-1))X1.0DO)/K)XQ 
END IF 

CUMSUM = CUMSUM + PTEMP1 
CONTINUE 

CUMSUM = CUMSUM X PP 
P(I,J) = L OGC CUMSUM) 

PRINT X, »I=',I, »J=‘, J, *LN= , ,P(I,J), *CUMSUM= ', CUMSUM 
CONTINUE 
CONTINUE 
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225 

200 



HI L AMD = -DINF 
DO 200 1=1, NCOMP 

DO 225 J = 1 , SPBNDC I ) 

HLTEMP = ( P(I,J) - P(I,0) ) / 
IF (HLTEMP .GT. HILAMD) HILAMD 
CONTINUE 
CONTINUE 

PRINT X, 1 HI LAMB DA = 1 , HILAMD 

RETURN 

END 



WCI, J) 

= HLTEMP 



SUBROUTINE BOUNDC DI NF, P , W, NCOMP , MAXSPR, SPBND, MAXWT , HI L AMD, NSPARE, 
XZL0W,ZUP, WT, BSPARE) 

INTEGER NCOMP, BSPARE( NCOMP), MAXSPR, SPBNDC NCOMP ),NSPARE( NCOMP) 

REALX8 BSTLAM, ZLOW,WT,P( NCOMP, 0: MAXSPR) ,H( NCOMP, 0 : MAXSPR) , DINF, 

XMAXWT, HILAMD, EPS, LFTEND, RTEND, ZLBEST, LAMBDA, ZUP 



EPS = HILAMD/1 . ODOA 
LFTEND = O.ODO 
' RTEND = 1 . 1D0XHILAMD 
ZLBEST = -DINF 

50 LAMBDA = (LFTEND + RTEND)/2.0D0 

CALL MAXFCT( DINF, P,W, NCOMP, MAXSPR, SPBND, MAXWT, LAMBDA, NSPARE, ZLOW, 
*ZUP,WT) 

IF (WT .LE. MAXWT) THEN 
RTEND = LAMBDA 
IF (ZLOW .GE. ZLBEST) THEN 
ZLBEST = ZLOW 
BSTLAM = LAM3DA 
END IF 
ELSE 

LFTEND = LAMBDA 
END IF 

IF ( (RTEND - LFTEND) . GT . EPS ) GOTO 50 

CALL MAXFEAC DINF, P,W, NCOMP, MAXSPR, SPBND, BSTLAM, BSPARE, WT, ZLOW) 

RETURN 

END 
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SUBROUTINE MAXFCTC DINF, P , W, NCOMP , MAXSPR, SPBND, MAXWT , LAMBDA , NSPARE, 
*ZLOW,ZUP,WT) 

INTEGER NCOMP, NSPAREC NCOMP) , MAXSPR, SPBNDC NCOMP ) , INDEX 

REALX8 L AMBDA,ZLOW,WT, PC NCOMP, 0 : MAXSPR ) ,W( NCOMP, 0 : MAXSPR), DINF, 

XCTOTAL , CMPMAX, ZVAL J , OBJFCT ,OBJVAL,ZUP, MAXWT 



WT = O.ODO 
CTOTAl = O.ODO 
ZLOW = O.ODO 

DO 25 1=1, NCOMP 
CMPMAX = -DINF 
INDEX = 0 

DO 50 J=0,SPBND(I) 

ZVAL J = P(I,J) - LAMBDAXWCI, J) 

03JVAL = P(I,J) 

IF ( ZVAL J .GE. CMPMAX) THEN 
NSPARE(I) = J 
CMPMAX = ZVAL J 
OBJFCT = OBJVAL 
INDEX = J 
END IF 

50 CONTINUE 

CTOTAL = CTOTAL + CMPMAX 
ZLOW = ZLOW + OBJFCT 
WT = WT + WC I, INDEX) 

25 CONTINUE 

ZUP = CTOTAL + LAMBDAXMAXWT 

RETURN 

END 

SUBROUTINE MAXFEAC DINF, P, W, NCOMP , MAXSPR, SPBND, BSTLAM, BSPARE , WT , 
XZLOW) 

INTEGER NCOMP, BSPAREC NCOMP ) , MAXSPR, SPBNDC NCOMP) , INDEX 

REAL 5(8 BSTLAM, ZLOW, WT, PC NCOMP, 0 : MAXSPR ),W( NCOMP, 0 : MAXSPR), DINF, 
XCMPMAX, ZVAL J , OBJFCT, OBJVAL 



WT = O.ODO 
ZLOW = O.ODO 

DO 25 1=1, NCOMP 
CMPMAX = -DINF 
INDEX = 0 

DO 50 J = 0 , SPBNDC I ) 

ZVAL J = P(I,J) - BSTL AM*W( I , J ) 
OBJVAL = P(I,J) 

IF C ZVAL J .GE. CMPMAX) THEN 
BSPAREC I ) = J 
CMPMAX = ZVAL J 
OBJFCT = OBJVAL 
INDEX = J 
END IF 

50 CONTINUE 

ZLOW = ZLOW + OBJFCT 
WT = WT + WC I, INDEX) 

25 CONTINUE 

RETURN 

END 
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SUBROUTINE HERIST ( DINF, P , W, NCOMP , MAXSPR, ZHLOW, WT , MAXWT , SPBND, 
XBSPARE, LARRAY) 

INTEGER NCOMP, B SPARE (NCOMP ) , MAXSPR, SPBND (NCOMP ), INDEX 

REALMS DINF , P ( NCOMP , 0 : MAXSPR) , W( NCOMP , 0 : MAXSPR ) , WT1 , WT , MAXWT , 
XXSSWT, CHKVAL , LARRAYCNCOMP ) , NUMER, DENOM, LBEST , ZHLOW 

CHKVAL = DINF/1.1D0 
WT 1 = WT 

25 XSSWT = MAXWT - WT1 
DO 50 1=1, NCOMP 

IF (CBSPARE(I) .GE. SPBNDCI)) .OR. (PC I , (BSPAREC I ) + l ) ) .GT. 
XCHKVAL)) THEN 

LARRAY ( I ) = -DINF 
GOTO 50 
END IF 

NUMER = P(I,(BSPARE(I)+1)) - P ( I , BSPAREC I ) ) 

DENOM = WCI, CBSPARECD + l) ) - W( I , BSPAREC I ) ) ■ 

IF ((DENOM .GT. XSSWT) .OR. (NUMER .LE. 1.0D-10)) THEN 
LARRAY ( I ) = -DINF 
GOTO 50 
ELSE 

LARRAY ( I ) = NUMER/DENOM 
END IF 

IF ( LARRAY ( I ) .LT. 0.0D0) LARRAY ( I ) = -LARRAY ( I ) 

50 CONTINUE 

LBEST = -DINF 
DO 75 1=1, NCOMP 

IF ( LARRAY ( I ) .GT. LBEST) THEN 
LBEST = LARRAY ( I ) 

INDEX = I 
END IF 

75 CONTINUE 

IF (LBEST .LT. (-CHKVAL)) GOTO 100 

WT1 = WT1 - W( INDEX, BSPAREC INDEX)) + W( INDEX, ( BSPAREC INDEX)+1 ) ) 
BSPARE( INDEX) = BSPAREC INDEX) + 1 
IF ( WT1 .LE. MAXWT) GOTO 25 
100 ZHLOW = 0 . 0D0 
WT = 0.0D0 
DO 125 1=1, NCOMP 

ZHLOW = ZHLOW + P( I , BSPAREC I) ) 

WT = WT + WCI, BSPAREC I ) ) 

125 CONTINUE 

RETURN 

END 
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APPENDIX B 



SAMPLE INPUT FILES 



PRINCIPAL END ITEM FILE 



M60A1 5 A 60.0 
LVTP7A1 12 11 40.0 
M101A1 5 4 2000.0 



TABLE B-l 

PRINCIPAL END ITEM FILE SPECIFICATIONS 



LINE 


COLUMN 


DATA DESCRIPTION 


DATA TYPE 


All 


1-7 


Nomenclature of end item 


Character 




8-11 


Blank 






12-14 


Quantity of end item 


Integer 




15-18 


Blank 






19-21 


Number of end items 
required operational 


Integer 




22-25 


Blank 






26-31 


Mission duration 


Real*8 
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SECONDARY REPARABLE COMPONENT FILE 



28 15001245387 


1 


2091.0 


7 


2000 . 0 


5 


0.34000 


58.6903 


2520002241867 


1 


2091 . 0 


8 


1000 . 0 


5 


1 . 00000 


1 . 0000 


2815002395819 


1 


2091 . 0 


3 


500 . 0 


5 


1 . 00000 


1 . 0000 


2520001549652 


1 


2091 . 0 


5 


500.0 


5 


0 .10110 


13.1850 


2530000886657 


1 


2091.0 


1 


100 . 0 


10 


1 .00000 


1 .0000 


2520001184942 


1 


2091 . 0 


1 


250.0 


5 


1 . 00000 


1 . 0000 


2920010135802 


1 


2091.0 


1 


500.0 


10 


1 .00000 


1 . 0000 


1010010703303 


1 


2091.0 


1 


1000.0 


10 


1 . 00000 


1 . 0000 


1025019388105 


1 


2091.0 


5 


100.0 


5 


1 . 00000 


1 . 0000 


1010010720397 


1 


2091.0 


3 


2000.0 


5 


1 . 00000 


1 . 0000 


1025010708993 


1 


2091.0 


1 


1000.0 


5 


1 .00000 


1 . 0000 


2815004303480 


2 


3253.0 


8 


500 . 0 


5 


0 . 44330 


142.9404 


2520003973384 


2 


3258.0 


23 


1000.0 


5 


0.69700 


69 . 0598 


2520001443385 


2 


3258.0 


14 


2000.0 


5 


1 . 00000 


1 . 0000 


5805014593410 


2 


3258.0 


10 


250.0 


5 


0.20310 


27 .7348 


2520008949532 


2 


3258.0 


3 


500.0 


5 


1.00000 


1 . 0000 


2910011714636 


2 


3258.0 


1 


100.0 


10 


1.00000 


1 .0000 


2530011509757 


2 


3258.0 


4 


1000.0 


10 


1 . 00000 


1 . 0000 


5865012207848 


2 


3253.0 


2 


100 . 0 


10 


1 . 00000 


1 . 0000 


4320012035660 


2 


3258.0 


3 


500 . 0 


10 


1 . 00000 


1 . 0000 


2920002317276 


2 


3258.0 


1 


500 . 0 


10 


1 . 00000 


1 . 0000 


2530000386650 


2 


3258.0 


1 


250.0 


5 


1 . 00000 


1 . 0000 


1005011355059 


5 


20996 . 0 


3 


1000.0 


5 


1.00000 


1 . 0000 


1005011086434 


3 


20996 . 0 


11 


250.0 


5 


1 . 00000 


1 . 0000 


1005011457709 


3 


20996.0 


7 


500.0 


5 


1.00000 


1 . 0000 
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SECONDARY REPARABLE/COMPONENT FILE SPECIFICATIONS 



LINE 


COLUMN 


DATA DESCRIPTION 


DATA TYPE 


All 


1-13 


National Stock Number ( NSN ) 


Character 




14-15 


Blank 






16-17 


Identification Number 
(1 for M60A1, 2 for LVTP7A1 , 
3 for M101A1 ) 


Integer 




18-19 


Blank 






20-28 


Total exposure time for 
component 


Real*8 




29-30 


Blank 






31-34 


Number of observed failures 
for component 


Integer 




35-36 


Blank 






37-42 


Component shipping weight 


Real*8 




43-44 


Blank 






45-47 


Maximum possible stockage 
for component 


Integer 




48-49 


Blank 






50-56 


Shape parameter of prior 
distribution for component 
failure rate 


Real*8 




57-58 


Blank 






59-66 


Scale parameter of prior 
distribution for component 
failure rate 


Real*8 
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APPENDIX C 



SAMPLE OUTPUT FILE 

OBJECTIVE FUNCTION VALUE IS s -0 . 8S20263153D+00 

THIS SOLUTION IS WITHIN 0 . 476 1377367D+00 % OF THE OPTIMAL 



PCMIS5I0N SUCCESS) 


= 0.41394 




NSN 


END ITEM 


SPARE 1 


2815001245387 


M60A1 


2 


2520002241867 


M60A1 


3 


2S 150 02395819 


M60A1 


2 


2520001549632 


M60A1 


2 


2530000886657 


M6 0A1 


3 


2520001184942 


M60A1 


2 


2920010135802 


M6 0A1 


1 


1010010703803 


M60A1 


1 


1025019388105 


M60A1 


5 


1C10010720397 


M60A1 


1 


1025010708993 


M60A1 


1 


2815004303480 


LVTP7A1 


4 


2520003973334 


LVTP7A1 


5 


2520001443335 


LVTP7A1 


4 


5305014593410 


LVTP7A1 


5 


2520008949532 


LVTP7A1 


2 


2910011714636 


LVTP7A1 


3 


2530011509757 


LVTP7A1 


2 


5365012207348 


LVTP7A1 


3 


4320012035660 


LVTP7A1 


o 


2920002317276 


LVTP7A1 


1 


2530000886650 


LVTP7A1 


2 


1005011855059 


M101A1 


4 


1005011086434 


M101A1 


5 


1005011457709 


M101A1 


5 
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GAMS SOURCE CODE 



SETS 

I end item /E1875, E03A6 , E06A0/ 

J component /1*25/ 

K spare stockage level /50*S10/ 

L dummy index /D0*D11/; 

SCALARS 

WTCAP total weight capacity (lbs) /AAAOO/ 
PSUC probability of mission completion; 

PARAMETERS 

EQDENS(I) density of end item i 
/E1S75 5 

E0SA6 12 

EQoAO 5/; 



PARAMETERS 

NOPER(I) number end item i required operational during mission 
/E1375 A 

EOS A6 11 

E06A0 A/; 



PARAMETERS 

DURMIS(I) 
/ El 87 5 
E08A6 
E06 AO 



mission duration for end item i (hours rounds miles) 
60 
AO 

2000 /; 



PARAMETERS 

SHIPWT ( J ) 
/I 
2 
3 
A 

5 

6 

7 

8 

9 

10 
11 
12 
13 
1A 

15 

16 

17 

18 

19 

20 
21 
22 
23 
2A 
25 



shipping weight of component 
2000 
1000 
500 
500 
100 
250 
500 
1000 
100 
2000 
1000 
500 
1000 
2000 
250 
500 
100 
1000 
100 
500 
500 
250 
1000 
250 
500/; 



(lbs) 
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PARAMETERS 

SPRMAXC J ) 
/I 
2 
3 
A 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
25 

24 

25 

PARAMETERS 
N0BS( J ) 

/I 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 



max number of spares allowed for component 

5 

5 

5 

5 

10 

5 

10 

10 

5 

5 

5 

5 

5 

5 

5 

5 

10 

10 

10 

10 

10 

5 

5 

5 

5 /; 



number of data points observed for component 

7 

8 
3 
5 
1 
1 
1 
1 
5 
3 
1 
8 

23 

14 

10 

3 
1 

4 
2 
3 
1 
1 
3 

11 

7 /; 



J 



3 
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PARAMETERS 
ALPHA( J ) 
/I 
2 
3 
A 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 



shape paramtr for Gamma prior distribution 
0.34 
1 
1 

0.1011 

1 

1 

1 

1 

1 

1 

1 

0.4433 

0.6970 

1 

0.2031 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 /; 



PARAMETERS 
BETA( J ) 
/I 
2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 



scale paramtr for Gamma prior distribution 
58.6903 
1 
1 

13.1850 

1 

1 

1 

i 

1 

1 

1 

142.9404 

69.0598 

1 

27.7348 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 /; 
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PARAMETERS 

0K(I,J) equals 1 if component j is a component of end item 1 
/E1875.C1X11) 1 

£0846.(12X22) 1 

E0640 . (23X25) 1/; 

PARAMETERS 

TXPOSE(I) total observed exposure for end item i ( hrs rds mi) 
/E1S75 2091 

E0846 3258 

E06 40 20996/; 

PARAMETER 

LMT(I,J) max number of failures for component j; 

LMT(I,J)$( 0K(I,J) EQ 1 ) = EQDEMS(I) - NOPERd); 



SETS 

OKLVL 1 ( I > J , K) feasible spare stockage levels for component j; 
0KLVL1(I,J,K)$(( 0RD(K)-1 LE SPRMAX(J)) AND (0K(I,J) EQ D) = YES; 

X The following set is necessary to obtain the cumulative distribution 
x function of the negative binomial density. 

SETS 

0KLVL2(I, J,K,l); 

0KLVL2(I, J,K,L)$( (ORD(K)-l LE SPRMAX(J)) AND (0K(I,J)EQ 1) AND 
(ORD(L)-l LE LMT(I,J)+0RD(K)-1) ) = YES; 



x 

x 



DISPLAY 0KLVL2; 

OPTIONS DECIMALS = 8; 

PARAMETER 

P(I,J) negative binomial paramtr 

Q(I,J) one minus the p paramtr 

PP(I»J) p paramtr raised to power of nobs plus aloha 

TWT(I,J»K) total snipping wt of k spares of component j 

C(I,J,K) cumulative density for negative binomial 

F(I,J,K) natural logarithm of cumulative density 

PTMP( I , J , K/ L ) temporary paramtr needed to sum cdf recursively; 

P( I » J )$( 0K( I / J ) EQ 1) = 

(TXPOSE(I) + BETA(J))/((EQDENS(I)XDURMIS(I)) + TXPOSE(I) + BETA(J)); 
Q(I,J)$(0K(I,J) EQ 1) = 1 - P ( I / J ) ; 

PP(I,J)$(0K(I,J) EQ 1) = P(I,J) xx ( NOBS ( J ) + ALPHA(J)); 

TWT(I, J,K)$0KLVL1(I,J,K) = SHIPWT(J) x (ORD(K) - 1); 

PTMP(I, J,K,"DO n )$OKLVLl(I, J,K) = 1; 

The following LOOP function is used to calculate the cdf from the 
negative binomial pdf recursively. 
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LOOP ( L , PTMPCI, J,K,L+1)$0KLVL2(I,J,K,L) = 

QCI, J)X( (NOBS(J)+ALPHA(J)+ORD(L)-l)/(ORD(L)))*PTMP(I, J,K,L) ); 

CCI, J,K)$0KLVL1CI, J,K) = 

SUM(L$((ORD(L)-l) LE ( LMTCI, J )+ORD(K)-l) ) , PTMPC I, J , K, L ) ) *PP( I . J ) ; 

FCI, J,K)$0KLVL1CI, J,K) = LOG( C(I,J,K) ); 

/ 

DISPLAY F; 

VARIABLES 

X(I,J,K) 1 if stock k spares for component j of end item i 
MOE natural log of probability of mission completion; 

BINARY VARIABLE X; 

EQUATIONS 

WTCONST weight constraint equation 

SELECT(I,J) select exactly 1 of the possible stockage levels 
MOEDEF measure of effectiveness definition; 

WTCONST.. SUMCCI, J , K)$0KLVL1 C I , J,K), TWT(I,J,K)XX(I,J,K)) =L= WTCAP; 

SELECT ( I , J ) $ ( 0K(I,J) EQ 1 ).. 

SUM( K$ OKLVL 1 ( I , J , K) , X(I,J,K) ) =E= 1; 

MOEDEF.. SUM((I,J,K)$0KLVL1(I,J,K), F( I , J , K ) 3€XC I , J , K) ) =E=M0E; 

MODEL STOCK /ALL/; 

OPTIONS OPTCR = 0.0015, OPTCA = 0; 

SOLVE STOCK USING MIP MAXIMIZING MOE; 

PSUC = EXP ( MOE.L ) 

DISPLAY X.L, MOE.L, PSUC; 
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